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The present paper extends the modern theory of tunneling transport to finite temperatures and 
presents an application for molecular devices made of alkyl chains connected to silicon nano- wires, 
mapping their transport characteristics as functions of temperature and alkyl chain's length. Based 
on these calculations and on the analytic theory, it is found that the tunneling decay constant is 
determined not by the Fermi level, but by the edge of the valence or conductance band, whichever is 
closer to the Fermi level. Further insight is provided by mapping the evanescent transport channels of 
the alkyl chains and few other physical quantities appearing in the analytic formula for conductance. 
A good qualitative agreement with the experimental data is obtained. 
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Our understanding of charge transport through single 
molecules is growing rapidly due to advances on both ex- 
perimental and theoretical fronts. This is especially true 
for molecular devices connected to metallic leads. How- 
ever, the theoretical understanding of transport in de- 
vices connected to semi-conducting leads is happening at 
a much slower pace. The number of atomistic simulations 
for this class of devices is quite limited. 1 6 This is some- 
what surprising giving the fact that today's electronic 
devices are mostly made of semi-condicting materials. It 
is also surprising because the self-assembled alkyl chains 
on silicon surfaces were among the first molecular elec- 
tronic devices to display highly reproducible transport 
and characterization measurements. 

The past few years have seen the development of what 
we call the modern theory of tunneling transport 
This theory connects the tunneling exponent /3 to the 
complex band structure of the insulating chains, an ap- 
proach that goes well beyond the effective mass treat- 
ments. In Refs . fTTI and \12\ one of the authors contributed 
to the theory by deriving an explicit expression for the 
contact conductance. The theory has been applied ex- 
tensively to devices with metallic leads, but never to de- 
vices with semi-conducting leads. For the later case, the 
following questions are open: a) At what energy should 
be j3 evaluated? b) Which states determine the contact 
conductance? c) How does the temperature enter in all 
of these? The goal of this paper is to clarified all these 
issues and to show the theory at work for the electronic 
devices discussed below. 

Electronic devices made of alkyl chains self-assembled 
directly on silicon surfaces present several key features 
allowing accurate and reproducible transport measure- 
ments and characterizations JI^Md p or example, silicon 
present s a sm ooth, reproducible and well-controlled solid 
surfacePH221 Also, the well-prepared Si-C bonded alkyl 
monolayers are stable, 23 and suitably long alkyl chains 
form well-organized, nearly pinhole-free monolayers 
The highly covalent Si-C bond, unlike the metal-thiolate 
or metal-amine bond etc., is a natural continuation of the 
alkyl chain, a feature that minimizes potential disconti- 
nuities at the contacts! 13 * 18 * 19 * The absence of a linker 



between Si and the alkyl chain allows the direct study 
of the electrical properties of the monolayer itself. The 
systems themselves are interesting because the choice of 
Si electrodes makes the devices relevant to Si-based mi- 
croelectronics. 

The electronic structure of self- assembled alkyl chains 
on Si(lll) surface was studied theoretically in Ref. [26] 
No such study exists for the charge transport properties 
of these devices. Using geometries t aken from experi- 
ments and previous theoretical workspE! Ref. [26] pro- 
vides a good correlation between the theoretically com- 
puted and experimentally measured ultraviolet photo- 
electron spectra (UPS) and inverse photoemission spec- 
troscopy (IPES). To obtain simultaneous agreement with 
the two sets of data, the theoretical Si insulating gap had 
to be properly adjusted by a rigid shift of the conduction 
states and the valence states had to be slightly stretched. 
These corrections were mainly necessary to correct for 
the smaller theoretical Si gap, a well known shortcoming 
of the current Density Functional Theory (DFT) kernels. 
However, the electronic structure around the band edges 
seemed to be reasonably well reproduced by these cal- 
culations, leading us to believe that DFT is appropriate 
for studying the transport when the Fermi level is close 
to the valence or conduction bands. Another interesting 
output of the calculations of Ref. [26] was that the lateral 
dispersion of the energy bands of the alkyls was weak, 
implying that the electronic states are localized on the 
chains even at relatively high monolayers coverage. This 
suggests that the transport properties of the monoloyers 
can be analyzed using just a single chain. 

The experiments mentioned above involve long alkyl 
chains, typically containing more than 12 methyl groups. 
We compute the transport characteristics of devices with 
only up to 8 methyl groups, but we show that the systems 
already reached the asymptotic regime and therefore we 
can extrapolate the conductance to larger devices. To 
compute the tunneling conductance of such large molec- 
ular devices, we use the asymptotic expression of the 
tunneling conductance derived in Refs. [11] and [12] We 
extend these derivations to finite temperature and then 
apply the new formalism to study the tunneling proper- 
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FIG. 1: (Color online) Illustration of a typical device con- 
sidered in this paper. It contains an alkyl chain covalently 
bonded to Si nano-wire leads. The figure indicates the unit 
cell that is repeated to obtain the periodic potential Vq. It 
also defines the length L of the chain. 



ties of single alkyl chains connected to semi-conducting 
Si nano- wires. The Si nano-wires were chosen because 
of certain computational advantages. Due to the lat- 
eral confinement of the electrons in Si nano-wires, the 
electronic structure of the electrodes quantitatively and 
qualitatively differs from that of bulk Si. Nevertheless, 
the present simulations allow us to map and rationalize 
the dependence of the transport characteristics on tem- 
perature and alkyl chain's length and to make qualitative 
comparisons with the experiments. 



I. THEORETICAL FRAMEWORK 

We consider a charge transport experiment involving 
a device made of an insulating molecular chain attached 
to semi-conducting leads (see Fig. [I]). The distinction 
between insulating and semi-conducting chains is that 
the former have much larger energy gaps. The device is 
assumed oriented along the z axis. The charge current is 
driven by a small time oscillating electric field E ext (r)e 2a;t , 
whose effects are treated in the linear response regime. 
The dc regime is obtained by letting the frequency of the 
oscillations go to zero. The existence of a steady state is 
implicitly assumed. 

The device will not conduct electricity at zero temper- 
ature. However, we can simplify the discussion if we start 
from the zero temperature formalism and assume metal- 
lic leads. For this case it was recently established thatp^ 
within the T ime Dependent Current-Density Functional 
Theoryp^l ^he linear conductance of a molecular device 
is given by: 
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Here, f is the matrix kernel 
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E^ n (r) being the dynamical response of the electrons to 
the external perturbation,^^ and <jj|(r, r'\ e F ) is the 
local conductivity tensor of the non-interacting Kohn- 
Sham electrons. r_\_ and r' ± denote the coordinates of 
two arbitrary planes normal to the axis of the device. 
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FIG. 2: (Color online) The atomic configuration of the Si 
wires. The left diagram shows a view along the wire and the 
next two show views from the sides. 



In the tunneling experiments, the currents are ex- 
tremely small and the non-adiabatic effects are expected 
to be negligible. In this case, we can neglect the term 
a KS *J^ and the zero temperature linear conductance 



takes the simple form 



G = J dr_L J dr' ± a^(r,r f ;e ¥ 



(3) 



This expression was previously derived in Refs. I3T1 and l32l 
for non-interacting electrons and it is equivalent to the 
Landauer formula applied to the non-interacting Kohn- 
Sham electrons. Therefore, Refs . ITTI and [T2l gives a formal 
justification why Landauer formalism is appropriate for 
treating the off-resonance tunneling transport of inter- 
acting electrons. 

Refs. I3l1 and l32l also established the finite temperature 
expression for the conductance of non-interacting elec- 
trons. We arrived at the same expression after including 
the finite temperature in the derivations of Refs. [TTland 
[12] that lead to Eq. [3j This expression is: 

G= 2 S I de[-fUe ~ e F )} T(e) (4) 



with 
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Here, /p D (e) is the derivative of the Fermi-Dirac distri- 
bution: J FD (e) = (1 + e^ 6 ) -1 [l//3=Boltzmann constant 
x Temperature]. The above expression applies equally 
well to devices with metallic or semi-conducting leads. 
To avoid confusion, we use r to denote the temperature. 

The exact asymptotic expression, in the limit of long 
alkyl chains, of T(e) was computed in Refs. [11] and [12] 
Let us give a brief discussion of this result. Assume that 
a self-consistent Kohn-Sham calculation has been com- 
pleted for the entire device. As we shall see later in this 
paper, to a high degree of accuracy, the effective poten- 
tial V eH of the entire molecular device can be decomposed 
into a perfectly periodic piece Vb, extending from nega- 
tive to positive infinity, and a difference AV=V eii — Vq. 
The periodic potential Vq is constructed by periodically 
repeating the effective potential between —b/2 and 6/2 in 
the middle of the chain, as indicated in Fig. [I] It is useful 
to regard the self-consistent Kohn-Sham Hamiltonian of 
the entire device as a periodic Hamiltonian, 
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FIG. 3: The band structure (left) and density of states (right) 
of the infinite Si nano-wire. The energy units are eV. 



strongly perturbed by the leads via the potential AV, 
which is localized on the electrodes and is practically zero 
inside the alkyl chain. The effective Hamiltonian of the 
entire system is then 



F = # + AV L (r)+AV R (r), 



(7) 



where we divided AV into the left and right parts relative 
to the mid point of the device. 

The Bloch functions for the periodic potential Vo are 
usually characterized by a /c-vector and a band index. 
Since here we consider energies inside the insulating gap 
of Vo, the band index losses its precise definition and it is 
more useful to drop this index and see the k- vector of the 
Bloch functions as defined on a Riemann surface. This 
Riemann surface was described extensively in Ref. |33j It 
turns out that the electron band dispersion e& can be 
extended to the whole Riemann surface and if we pick 
an arbitrary energy (possibly complex like e F + iS) then 
the equation e& = e has a sequence {k a } of infinite solu- 
tions. We are interested in the case when e = e F , which is 
located inside the insulating gap of the alkyl chain. For 
this reason, all the fc a 's have finite imaginary parts, which 
are always taken with positive sign. For each such fc a , 
there is an evanescent Bloch solution ipk a (r) exponen- 
tially decaying to zero when z goes to positive infinity, 
and an evanescent Bloch solution ^_/ Ca (r) exponentially 
decaying to zero when z goes to negative infinity. The 
exponential decay factor of the evanescent solutions is 
given by Im[fc a ]. The alkyl chains conduct through all 
the evanescent channels, but for long chains the charge 
transport mainly occurs through the channel with min- 
imum Im[fc a ], the contributions from the other channels 
being exponentially small. Therefore, we can focus on 
the channel with the lowest Im[fc a ] and for this reason 
we drop the index a from now on. 

To compute the transmission of the whole device, writ- 
ten in Eq. [8j one can place the normal planes r± and 
in the middle of the device, in which case one only needs 



to understand how the evanescent channels reflect at the 
contacts on to the potentials V^ /R (r). This has been ac- 
complished in Refs. [TT1 and [T2l by using a Green function 
formalism. The asymptotic expression of T(e), in the 
limit of long alkyl chains, was found to be: 



T(e) = e L e R e- 2/m M L . 

The coefficients are given by: 
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with r and r' measured from the left end of the chain. 
Similarly 
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with r and r' measured from the right end of the chain. 
In both expressions, p eF (r,r') is the spectral operator: 
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The sum goes over all Kohn-Sham orbitals 4> e (r). The 
symbol W((j), ip) stands for the generalize Wronskian: 



W((j)^) = / dr± 4>(r)d z il)(r 



(12) 



The expressions for coefficients are independent of how 
one normalizes the evanescent Bloch functions. 

A closer look at the integrands of Eqs . |9| and [l0| reveals 
that ^_fc(r)AV^(r) and i/jk(r)AV K (r) are exponentially 
localized near the contacts. Therefore, to evaluate the 
6 coefficients, one only needs a converged spectral op- 
erator near the contacts, which can be obtained from a 
standard electronic structure calculation which includes 
long enough leads. 

The asymptotic expression of G becomes virtually ex- 
act for alky Is containing more than 4 methyl groups.^ 
The simplification brought in by this expression allows us 
to compute the tunneling conductances of the extremely 
large molecular devices presented in this paper without 
truncating the Hilbert space. We use the same number of 
basis functions (better said grid points) in the transport 
calculation as in the self-consistent calculation, which for 
one device is larger than 2.1 xlO 6 . Most of the transport 
calculations completed so far use a number of basis func- 
tions that is orders of magnitude smaller than the ones 
used in this paper. 

Besides the computational advantages, the analytic ex- 
pression for G allows us to get an unprecedented in- 
side look into the transport characteristics of the devices. 
Since the expressions for the coefficients are basically 
overlap integrals, each quantity inside the integrands re- 
veals one specific aspect. The evanescent Bloch functions 
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FIG. 4: (Color online) (a) The atomic configuration of the 
device, b) An iso-surface of V e a. c) The xy average of V e ff. d) 
An iso-surface of AV. e) The xy average of AV . The energy 
units in all panels are in eV. 



allows us to image the path that is important for tunnel- 
ing transport. A plot of '0_/ c (r)A 1 l/ L (r) and ^k(r)AV n (r) 
reveals the localization of the contact conductance, al- 
lowing us to determine which layers inside the leads are 
important for the tunneling transport. 

The formalism is applied to devices made of alkyl 
chains connected to Si nano- wires, as illustrated in Fig.[l] 
Various alkyl chains and Si nano-wires will be considered. 
All the broken Si bonds are passivated with hydrogen and 
all the bonds are kept in a perfect tetrahedral configu- 
ration. The Si atoms are arranged in the bulk configu- 
ration, the S-H bond length is fixed at 1.49 A and the 
Si-C bond at 1.75 A. For the alkyl chains, the C-C bond 
length is fixed at 1.54 A and the C-H bond length at 1.1 

A. 

The calculations are performed in a periodic supercell 
containing the alkyl chain and the Si nano-wires as illus- 
trated in Fig. [I] The equilibrium self-consistent Kohn- 
Sham calculations, the real and complex band structure 
calculations and the non-equilibrium transport calcula- 
tions were performed with a real space, pseudopoten- 
tial code based on finite differences. The same code was 
used for the calculations reported in Ref. |34] and [T2j We 
adopted a 5-point (in each direction) finite difference ap- 
proximation for the kinetic energy operator, and used 
a uniform space grid of spacing 0.19 A, sufficient for a 
good convergence of the electronic band structure. This 
grid is commensurate with the unit cell of the periodic 
alkyl chain, which is the reference system in our trans- 
port calculations. The largest device discussed in this 



paper requires the staggering number 81x81x329 of grid 
points and contains 1092 electrons. The lateral grid size 
is kept constant at 81x81 grid points or 15.4Axl5.4A. 

We adopted the Local Density Approximation (LDA) 
for exchange and correlation using the Perdew-Zunger 
(PZ) 35 interpolation of the numerical electron-gas data 
of Ceperley and Alder P*l We used Troullier-Martin norm- 
conserving pseudo-potentials 3 -^ for all the atomic species. 
The pseudopotentials for C and Si atoms had distinct s 
and p components and we took the p pseudo-potential as 
the local reference. Purely local pseudopotentials were 
used for the H atoms. 



II. THE SEMI-CONDUCTING ELECTRODES 

Although the present work was motivated by the trans- 
port measurements on self- assembled alkyls on Si(lll) 
surface, we cannot approach these systems at the mo- 
ment. Due to the tilting of the alkyl chains relative to the 
Si(lll) surface, a computation on the alkyl mononlayer 
will require a non-rectangular suppercell. To avoid such 
computational complication, we chose to represent the Si 
electrodes by Si nano-wires grown along the [110] direc- 
tion. The Si nano-wires are cut along the [111] surface 
where the alkyl chain is then attached. In principle, by 
taking the diameter of the wire larger and larger, we can 
model a low density alkyl monolayer. But on the other 
hand, since Si nano-wires have been fabricated for some 
time, the devices discussed in this work could very well 
be fabricated in the near future. For example, free stand- 
ing silicon wires grown along directions like [100], [110] 
or [112] have been already observed experiment allyP^ll 
Their electronic properties and their potential for novel 
electronic devices were discussed in Refs. I40H441 

In this work we consider the thinnest [110] hydrogen 
passivated Si nano-wires. The unit cell of the wire con- 
tains 16 Si and 12 H atoms. The atomic configuration 
of the Si nano- wire is illustrated in Fig. [2] from various 
view angles. The Si atoms are fixed at the bulk positions 
as previously discussed. This model structure is a fairly 
accurate representation of the relaxed configurations, the 
difference between the two being less than 1%. 42 

As already pointed out in Ref. ED the confinement of 
the electrons in thin Si wires leads to an electronic struc- 
ture that is quite different from the bulk Si. In Fig. [3] we 
show the band structure of our Si nano- wire and the cor- 
responding density of states. The band structure and the 
density of states near the band edges in Fig. [3] is similar 
those reported by previous calculation. 42 The nano-wire 
displays a direct gap of 1.75 eV, in line with the pre- 
vious theoretical works. 42 As one can see in Fig. [| the 
insulating gap is determined by highly dispersive bands. 
This bands contribute very little to the total density of 
states but, as we shall see, they have large effect on the 
transport. This particularity of the Si nano-wires makes 
the transport problem more interesting, but at the same 
time more difficult because it slows the convergence with 



FIG. 5: (Color online) The atomic configurations of the three molecular devices used to investigate the thermodynamic limit 
and the corresponding xy average of the local density of states. The energy unit is 1 eV. 
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FIG. 6: Real (right panel) and complex (left panel) band 
structures of the periodic potential Vb corresponding to the 
device shown in Fig. [4] Only the complex band with smallest 
Im[/c] is shown. The energy unit is 1 eV and the unit for Im[fc] 
is 1/Bohr. 



the leads size. 



III. ELECTRONIC STRUCTURE: GENERAL 
CONSIDERATIONS 

First, let us discuss the periodicity of the effective po- 
tential inside the alkyl chain, because this is the main 
assumption of our theory. The full effective potential has 
local and non-local parts. Since we work with norm- 
conserving pseudopotentials, the non-local part of the 
potential is not changed during the self-consistent calcu- 
lation. Therefore, the non-local part of the effective po- 
tential remains perfectly periodic inside the alkyl chain. 
In Fig. |4^b-c) we plot the local part of the effective po- 
tential for the entire molecular device, shown in Fig.[4ja). 
According to these graphs, the effective potential is pe- 
riodic inside the alkyl chain, starting from the very first 
methyl group. The periodic potential Vo is obtained by 



repeating the local and non-local parts of V eff inside the 
mid unit cell. If we subtract this Vb from V efi we ob- 
tain the local potential difference AV, which is plotted 
in Fig. [4jd-e). AV rapidly decays to zero away from the 
contacts and stays at zero for most of the alkyl chain's 
length. By construction, the non-local part of Ay is au- 
tomatically zero on the alkyl chain. We recall that AV 
appears explicitely in the expression of the © coefficients, 
therefore Fig. [4] is also important for understanding the 
results on the conductance. The above conclusions hold 
true for all the electronic devices considered in this work. 

The second issue we want to discuss is the fact that 
our molecular devices display a global insulating gap. 
This is evident in all the plots showing the local den- 
sity of states (see for example Fig. |5|. The gap of the Si 
nano-wire is clean of electronic states because the bro- 
ken Si bonds were carefully passivated with hydrogens. 
It is also important to notice that the gap is clean at the 
contacts, and this is again because we passivated the bro- 
ken Si bonds resulting from clipping the Si nano-wires. 
Referring to Fig. [5j one can distinctly notice the upper 
edge of the valence states for both Si and alkyl wires and 
the lower edge of the conduction states of the Si nano- 
wires, but only a faint mark of the conduction states 
edge of the alkyl chain can be observed. This is because 
we have computed only 75 un-occupied Kohn-Sham or- 
bitals. Since we do not have a good description of the 
conduction bands, in our transport calculations we as- 
sume p-doped wires and therefore a Fermi level that is 
close to the valence bands. 

The self-consistent calculations were performed at a 
finite temperature. In this case, we have a well defined 
Fermi level, but its position is irrelevant because it is 
highly sensitive to the presence of any impurity or defect, 
which will inherently occur in real devices. Because we 
include a fairly large vacuum region around the devices, 
the vacuum level is well defined and we can reference the 
electronic states from it. This is done throughout the 
paper. 

The real and complex band structures of the periodic 



6 



^ io~ 4 





r 1 rj. 


r\ 




T =0K 




I=300K 












r = 400 K 




T = 500 K 





-7 -6.5 -6 -5.5 -7 -6.5 -6 -5.5 



£ (eV) £ (eV) 

FIG. 7: (Color online) Plots of the finite temperature trans- 
mission computed at 0, 300, 400 and 500 K, for the three 
device shown in Fig. [5] The colors green, blue and red are 
used for devices (a), (b) and (c), respectively. The shaded 
region indicates the energy range of the valence states of the 
devices. The Fermi level of the device is located to the right 
of this region. 

potential Vb for the largest device of Fig. [5] are shown 
in Fig. [6j The complex bands remain almost unchanged 
when comparing one device to another. Apart from a 
rigid shift, both structures in Fig. [6] are similar to those 
reported in Ref. [45] for free standing alkyl chains. This 
tells us that, apart from a rigid shift coming from band 
alignment, Vo is similar to the effective potential of free 
standing chains. In Fig. [6] we show only the complex 
bands that are relevant to the tunneling. Note that the 
upper edge of the insulating gap practically terminates 
into the vacuum. This is a well known shortcoming of 
LDA, which places the bands too high in energy, though 
the shape of the bands is quite well reproduced, when 
compared to the experimentally resolved bands ™ For 
this reason, LDA gives a much smaller gap than the ex- 
perimentally measured one. Nevertheless, if one follows 
the complex band originating from the molecular valence 
band, 45 one will find that it connects to a band that is 
approximately 9 eV higher in energy, an energy difference 
that is close to the experimentally measured insulating 
gap of the alkyl chains. For this reason, the complex 
band shown in Fig. [6] is expected to have the right shape, 
thus leading to correct values of /?, provided the Fermi 
level is correctly predicted. 



IV. THERMODYNAMIC LIMIT 

In this section we discuss the convergence of our results 
with the size of the leads, which is always an important is- 
sue when working with finite systems) 4 ^ For our case, we 
studied three devices with finite Si nano-wire leads con- 




FIG. 8: (Color online) Plots of the finite temperature trans- 
mission computed at 300 K (blue), 400 K (gree) and 500 K 
(red) . The inset shows the values of the linear conductance at 
these three temperatures, with the Fermi level assumed to be 
position at the mark shown in the main graph. The contin- 

e F — e V 

uous line in the inset is a plot of the function ge ^ , with 
g adjusted till the best fit was obtained. The shaded region 
indicates the energy range of the valence states of the devices. 
The Fermi level of the device is located to the right of this 
region. 

taining a total of 9, 11 and 13 unit cells. These molecular 
devices are shown in Fig. [5] In the same figure, below 
the atomic structures, we show the corresponding xy av- 
erage local density of states: p av (z,e)=J p e (x,y, z)dxdy. 
The plots give a color map of pav(^e) in the plane of 
energy e (in eV) and of position along the device. For 
the density of states, we used a smearing factor S of 0.2 
eV. At the scale of these plots, the local density of states 
look completely converged, the plots showing no differ- 
ence when passing from one device to the other. It is also 
important to point out that the edges of the valence and 
conduction states for both Si nano-wire and alkyl chain 
do not change from one device to the next. 

Smearing factors of 0.1 to 0.2 eV are typical in trans- 
port calculations of T(e)PHere, however, since we per- 
form finite temperature calculations, we do not want the 
S smearing to interfere with the physical temperature 
smearing. Therefore, we must use a much smaller S than 
in the typical transport calculations. For example, the 
transport calculations at 300 K would require a 5 less 
than 0.026 eV and at 500 K a S less than 0.043 eV. On 
average, the energy level spacing in our devices is much 
less than 0.02 eV. Unfortunately, the transport is very 
sensitive to the electronic structure near the top of the 
valence band (assuming p-doped wires). This structure is 
determined by a very dispersive band therefore the spac- 
ing between the energy levels near the top of the valence 
band becomes much larger. For example, the average 
level spacing for the first five levels from the top of the 
valence band is 0.11 eV. 



FIG. 9: (Color online) The atomic configurations of the three molecular devices used to investigate the dependence on the 
number of monomers of the alkyl chain and the corresponding xy average of the local density of states. The energy unit is 1 
eV. 



The above criteria are only qualitative and we should 
look directly at the conductance results, which are shown 
in Fig. [7] for various temperatures. The K plot shows 
the transmission T(e) computed with a 5 smearing of 0.01 
eV, which was kept at this value throughout the paper. 
The next plots show the convolution of T(e) with the 
first derivative of the Fermi-Dirac distribution at different 
temperatures r: 

T(e,r) = |[-/; D ( e -e';r)]T(e')de'. (13) 

We call this quantity the finite temperature transmis- 
sion. The finite temperature conductance of the devices 
is then given by T(e,r)\ e=eF . T(e,r) can be accurately 
computed only in a limited energy window. The lower 
limit of this energy window is determined by the edge 
of the Vo potential's valence band, because our theory 
cannot be applied below this energy limit. The upper 
limit of the energy window is determined by the energy 
where the spurious Lorentian tails seen in the first panel 
of Fig. [7] become larger than the physical tails of the 
Fermi-Dirac distributions. Beyond this point, T(e,r) de- 
cays as a Lorentian, a spurious effect due to the finite S 
smearing used in the calculations. 

Looking at T(e, r) for different temperatures, we see 
the thermodynamic limit being accelerated as the tem- 
perature increases. The size dependence becomes quite 
small in the entire energy window, if one looks at the 
plots for 500 K. But the most important fact is that 
T(e,r) shows almost no size dependence in the energy 
window above the valence states, for all three tempera- 
tures. 



V. DEPENDENCE ON THE TEMPERATURE 

In this section we explore the behavior of the tunnel- 
ing conductance for different temperatures. In Fig. [8] we 



show T(e,r) at r=300, 400 and 500 K for the largest 
device shown in Fig. [5] For energies below the insulat- 
ing gap, the temperature effects would be negligible if we 
had a converged density of states. In our case, we do 
see a temperature effect below the insulating gap and, as 
expected, T(e,r) becomes smoother as the temperature 
increases. Inside the insulating gap, we see large varia- 
tions of T(e, r) with the temperature and these variations 
are physical. It is important to observe that the temper- 
ature dependence becomes weaker until it disappears as 
the Fermi level gets closer and closer to e v , the edge of 
the valence band. We will come back to this issue, after 
we map the dependence of conductance on the number 
of monomers, when we can put together the scenario de- 
scribing the experimental observations. 

To see how the conductance itself behaves with the 
temperature, we kept the Fermi level pinned at the same 
value as we varied the temperature. This value was cho- 
sen to be -5.55 eV, about 0.16 eV above the leads' va- 
lence band. From the main plot of Fig. [8j we collected 
the three values of the conductance corresponding to this 
e F , which are shown in the inset of the same figure. The 
continuous line in the inset of Fig. [8] represents the func- 

e F — e V 

tion ge ^ , which describes quite well the computed 
values for the conductance. This is precisely what the 
transport measurements for self-assembled alkyls on Si 
surface show. 

Few final notes for the Section are in place. The 
present calculations use the finite temperature DFT for- 
malism to account for temperature effects on the elec- 
tronic structure. However, since here we are discussing 
room temperatures (± 200 K), the electronic structure 
is practically unchanged from that at K. The issue 
of electron-phonon scattering inside the device could be 
raised. If the effect would have had any observable ef- 
fects, the behavior of G with the length of the insulating 
chain would deviate from the purely exponential tunnel- 
ing decay behavior. But we know from the extremely 
well controlled experiments on amine-linked alkyl chains 
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on gold leads'^ that the exponential tunneling decay be- 
havior is obeyed quite strictly [see also the discussion in 
Refs. l49l and [50] . There is no reason for things to be dif- 
ferent when the alkyl chains are attached to Si wires and 
therefore we conclude that the electron-phonon scatter- 
ing inside the device is negligible. It will definitely be- 
come an issue if one is concerned with the resistance of 
the long Si leads, which we are not. 

Another concerned could be structural molecular 
changes induced by temperature variations. The trans- 
port experiments P*"^ were never able to resolve any ef- 
fects coming from possible structural changes because, 
even if they occur, they are dwarfed by the large changes 
in the conductance already captured and discussed in 
Fig. [8} 
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FIG. 10: (Color online) Plots of the finite temperature trans- 
mission computed at 300 K and 500 K for the three devices 
shown in Fig. [9] The colors green, blue and red are used for 
devices (d), (e) and (f), respectively. The shaded region in- 
dicates the energy range of the valence states of the devices. 
The Fermi level of the device is located to the right of this 
region. 



VI. DEPENDENCE ON THE NUMBER OF 
MONOMERS 



For this part of our study, we work with the devices 
shown in Fig. [9j containing 4, 6 and 8 methyl groups. 
Similarly to Fig. |5j this figure illustrates the xy aver- 
aged local density of states for the three devices. The 
plots show a good alignment of the energy levels when 
comparing one device to the next. Also, the plots show 
almost no changes in local density of states near the con- 
tacts and inside the Si nano- wires, a fact that strongly 
suggests that the asymptotic regime was reached. 

We computed the finite temperature transmission for 
these devices at different temperatures and the results are 
shown in Fig. [Toj According to this data, if we fix our at- 
tention to one arbitrary energy above the valence band, 
the transmission T(e, r) decreases exponentially with the 
number N of monomers as T(e, r) oc e~@ N . Interestingly, 
/3 is practically independent of temperature. Since the 
expression for T(e, r) in Eq. [l3] is just the convolution 
of the zero temperature T(e') and the derivative of the 
Fermi Dirac distribution, T(e, r) is primarily determined 
by the e' that are close to the upper edge of the valence 
band. The contributions from the e' that are deeper into 
the valence band are exponentially suppressed. There- 
fore, we can not only understand the above observation, 
but we can actually tell what determines /3, something 
that eluded scientists for some time. For energies above 
the valence band, T(e, r) decays exponentially with the 
number of monomers and the exponential decay constant 
f3 is temperature independent and given by 2Im[/c]6 (see 
Eq. [8| with Im[k] evaluated at the top of the valence 
band. This is an important observation because it tells 
us that the tunneling decay constant in devices with semi- 
conducting leads is not determined by the Fermi level but 
rather by the edge of the valence band, for p-doped leads, 
and by the edge of the conduction band for n-doped leads. 
The present calculations give /3=0.73. 




FIG. 11: (Color online) Qualitative illustration of the local 
density of states for a real molecular device connected to a 
p-doped Si electrode. Typical energy values are also shown. 
The rectangle drawn with a dotted line represents the portion 
of the device that can we capture in our atomistic simulations. 



VII. BRIDGING WITH THE EXPERIMENT 

So far we have been careful not to place a firm Fermi 
level in our figures and we showed only results for the 
transmission rather conductance (with one exception), 
because determining the position of the Fermi level is 



beyond our atomistic simulations. In Fig. 11 we show a 
qualitative diagram of the local density of states in real 
devices, assuming p-doped Si electrodes. The box shown 
in this figure represents the part of the device that we 
can actually model with our atomistic simulations. In 
real devices, the positio n of t he Fermi level is determined 
by the impurity doping) 15 * 17 ! which is impossible to sim- 
ulate from first principles. Moreover, the doping triggers 
a band bending 18 that is also impossible to capture with 
any purely atomistic simulation. Therefore, e F — e v can- 
not be predicted by our calculations. We also want to 
mention that our transport simulations exclude the tun- 
neling through the Schottky barrier. However, this part 
is thoroughly understood and can be dealt with by well 
established theoretical treatments.^ 

In most measurements involving semi-conducting 
leads the experimental I-V curves display two dis- 
tinct regions corresponding to thermionic and tunneling 
regimes. For example in Ref. [I4j the thermionic regime 
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starts from zero bias potential and ends at about 0.6 
eV bias potential. It was argued that the main effect 
here comes from the flattening of the band bending as 
the bias is applied, bringing the alkyl chain's gap edge 
closer to the Fermi level (see Fig. 11). This part of the 
I-V curves was found to be sensitive to temperature and 
only slightly sensitive to the alkyl chain's length. Beyond 
the 0.6 V bias potential, the I-V curves display a linear 
behavior, with a slope that is highly dependent on the 
alkyl chain's length and almost independent of temper- 
ature. It was argued that bands are flat in this regime 
and that the I-V curves are determined by the tunneling 
through the alkyl chain. Almost as a rule, it was ob- 
served that the band bending in p-doped semiconductors 
is small and the tunneling regime starts from very low 
bias potentials.^ 



We argue that our finite temperature expression for the 
tunneling conductance can capture both regimes. We use 
our shorter alkyl chains to make qualitative contact with 
the experimental data mentioned above. To mimic the 
band flattening with the applied bias potential $, we fix 
the valence and conduction bands but consider a bias po- 
tential dependent Fermi level (in eV): e F (<I>)=e v +0.2— <I> 
for $<0.18 eV and e F ($)=e v +0.02 otherwise. This is 
consistent with the numbers shown in Fig. [TT] and as- 
sumes a band flattening that is linear with In reality, 
the band flattening will take place only on one side of 
the device, but we do not have the means to account for 
that fact at the moment. Only relatively small bias po- 
tentials will be considered so that we can apply the linear 
response theory on our molecular devices. In Fig. [12] we 
generate the I-V curves for the devices shown in Fig7|9j by 
using the formula /(<!>) =G(3>)-<!>, where the bias potential 
dependence in G(&) is induced by e F (<£). The expression 
for the current assumes that the potential <!> distributes 
entirely over the molecular device. For larger Schottky 
barriers we must take into the account that only a part 
of $ is actually distributed over the device. 



There is a strong qualitative agreement between our 
theoretical I-V curves and the experimental I-V curves 
for both n-doped and p-doped Si electrodes. The overall 
qualitative agreement is almost perfect with the experi- 
mental I-V curves for the p-doped electrodes of Ref. [T8l 
For the n-doped Si electrodes,^ the Schottky barrier is 
much higher and evidently our I-V curves are missing the 
effects of the this barrier. Nevertheless, the entire shape 
of the experimental I-V curves^ is well reproduced by 



our theoretical I-V curves. The insight of Fig. 12 shows 
two I-V curves computed at the two different tempera- 
tures considered in Ref. [I4j These theoretical I-V curves 
and the experimental I-V curves from Ref.[14j are in very 
good qualitative agreement. 




0.2 0.25 

<D(V) 

FIG. 12: (Color online) The predicted I-V curves for the de- 
vices (d), (e) and (f) computed at 330 K. The inset shows 
the predicted I-V curve computed at 330K and 255 K as in 
Ref. HI 



VIII. INSIGHT INTO THE TRANSPORT 
PROPERTIES 

Let us briefly describe how we computed the zero tem- 
perature conductance transmission. We limit the dis- 
cussion of this Section to the device (f) of Fig. [9] Given 
the particular complex band structure of the alkyl chains, 
the tunneling conductance is determined by just one com- 
plex band, the one with the smallest Im[fc]. This complex 
band is shown in Fig. [6] The complex band was obtained 
by varying continuously Im[fc] from to its maximum 
value, while keeping Re[fc]= 0. For each complex value 
of /c, the spectrum of the k dependent Hamiltonian: 



H k = -(V - ike z ) 2 + V + e 



-ik(z- 



-Xr,r'), (14) 



with periodic boundary conditions at z = ±6/2, was cal- 
culated and its eigenvalues ordered according to their real 
parts: Re[eife]<Re[e2fe]< •••• We focused, in particular, 
on the 6th and 7th eigenvalues e^k and ejk (which take 
real values, see Fig. [6| and their corresponding evanes- 
cent Bloch functions and ip^k- When lm[fc]=0, e§k 
and ejk coincide, respectively, with the top of the va- 
lence band and with the bottom of the conduction band 
of the periodic potential Vq. By increasing Im[fc], the 
two eigenvalues move towards each other until they meet 
and actually cross each other. At various values of Im[fc], 
we evaluated Eq. [8] for both e=€Qk and 6=67^, using the 
corresponding evanescent Bloch functions iJjQk and ip^k to 
compute the coefficients via formulas [9] and 10 The 
spectral kernel was computed directly from the Kohn- 
Sham orbit als of the full device as previously explained. 
The smearing coefficient S was fixed at 0.013 eV. 



The analytic expressions of Eqs. [8j [9] and 10 allows us 
to point several key aspects of the tunneling transport of 
our devices. The © coefficients are given by the overlap of 
the evanescent wavefunctions of the the alkyl chain, the 
potential perturbation AV and the spectral function of 
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FIG. 13: (Color online) Illustration of the overlap between the 
evanescent wavfunction V'fcfr) and (a) Ay and (b) p ev (r). 



the entire device. Therefore, we can gain valuable insight 
into the transport properties of our devices by mapping 
each of these three physical quantities and their overlap. 

A plot of the averaged local density of states (i.e. the 
diagonal part of the spectral operator) was already given 
in Figs. [5] and [9] Fig. [4] illustrates the decomposition 
of V eff in terms of Vo+AV. In Fig. [l3^a) we show the 
overlap between AV and the right evanescent function 
^fc(r) evaluated at e v . This overlap is further analyzed 
in Fig. [l4] Fig. 13 'b) shows the overlap of the evanescent 
wavefunction an d the local density of states p ev ( r ) 

at e v . The plot reveals a very strong overlap which ex- 
plains why the top valence states give such a large contri- 
bution to the transmission, while they contribute so little 
to the density of states (see Fig. |3|. We should comment 
here that the density of states at e v is determined by the 
highly dispersive band seen at and near the top of the 
valence band in Fig. [3j The character of the top valance 
bands in Si nano- wires have been thoroughly discussed in 
Ref. |4TJ and shown to originate from the lateral confine- 
ment of the electrons (the atomic character remains the 
same as in the bulk Si). In fact, the top valence bands 
are quite well reproduced by a jellium model confined in 
a cylinder of same radius as the wires. As the lateral size 
of the Si nano-wire is increased, the band structure con- 
verges to that of bulk Si, therefore the top valance bands 
become less dispersive and the associated wavefunctions 
less localized along the center of the nano-wire. Because 
of these two factors, that is, strong changes in dispersion 
and localization of the top valance bands, we expect a 
strong dependence of conductance on the lateral size of 
the Si nano-wire. 

We now return to the overlap between the evanescent 
wavefunctions ^±k{r) and AU L/R : 



* L/R (r) = ^(r)AK /R (r). 



(15) 



Since ^ T fe(r) decays exponentially inside the leads and 
AK/r decays extremely fast inside the alkyl chain, the 
two overlaps ^ L/R (r) are exponentially localized at the 
left/right contacts, respectively. As a consequence the 




(a) 





FIG. 14: (Color online) (a) Iso-surface plots (in green) of 
|^l/r(^)| imposed over the atomic structure of the device, 
(b) The xy average of |^ L / R (r)|. 



near the contacts, which is an important numerical ob- 
servation. But ^ L/R (r) also contains important physical 
information. An iso-surface plot of \I/ L/R (r) is shown in 
on top of the atomic structure of the device. 
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Fig._ 

This plot allows us to understand how different Si atoms 
contribute to the contact conductance. First thing we see 
in Fig. 
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'a) is that \P L / R (r) is very localized in the lat- 
eral direction therefore revealing a clear path inside the 
Si nano- wires which marks the region which is important 
to the contact conductance. Fig. [l4^a) also gives us a 
quantitative way to assess how deep into the Si nano- 
wires the overlap extends. For even better quantitative 
assessment, we show in Fig. 14 the xy average of ^ L/R (r). 



spectral operator in Eq. 11 is only needed in a region 



According to this graph, the contact conductance is de- 
termined by the first four Si layers of the nano- wires. 



IX. CONCLUSIONS 

In conclusion, we have extended the theory of the 
off-resonance tunneling transport to finite temperatures. 
The extension enables first-principles calculations of the 
charge transport in molecular devices connected to semi- 
conducting leads. We presented an application of the 
formalism to devices made of alkyl chains connected to 
Si nano-wires. 

The transport characteristics of the devices were char- 
acterized as functions of temperature and number of 
monomers in the alkyl chains. It was found that, for 
energy values above the valence states, the transmission 
decays exponentially with the number of monomers. The 
exponential decay constant is determined by the imagi- 
nary part of the complex wave- vector of the alkyl chains 
evaluated at the top of the valence band of the Si nano- 
wires. Using an empirical model of the band bending, we 
compared our results with experimental data and showed 
that the predicted I-V curves capture all the qualitative 
features of the experimental curves. 

The transport characteristics of the devices were fur- 
ther analyzed by mapping the physical quantities enter- 
ing in our analytic formula for transmission. We found a 
strong overlap between the evanescent waves of the alkyl 
chain and the local density of states at the very top of 
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the valence band, leading to an enhanced contribution of by the office of the Provost of Yeshiva University, 
these states to the transmission of the devices. Quanti- 
tative assessments of the localization of the contact con- 
ductance were also provided. 



Acknowledgments 



This research was supported by a Cottrell award from 
the Research Corporation for Science Advancement and 



1 T. Rakshit, G. C. Liang, A. W. Ghosh, and S. Datta, Nano 
Lett. 4, 1803 (2004). 

2 V. Meunier, W. Lu, J. Bernholc, and B. G. Sumpter, J. 
Phys.: Conf. Ser. 16, 283 (2005). 

3 W. Lu, V. Meunier, and J. Bernholc, Phys. Rev. Lett. 95, 
206805 (2005). 

4 S. Y. Quek, J. B. Neaton, M. S. Hybertsen, E. Kaxiras, 
and S. G. Louie, Phys. Rev. Lett. 98, 066807 (2007). 

5 G. Kirczenow, P. G. Piva, and R. A. Wolkow, Phys. Rev. 
B 72, 245306 (2005). 

6 G. Kirczenow, P. G. Piva, and R. A. Wolkow, Phys. Rev. 
B 80, 035309 (2009). 

7 P. Mavropoulos, N. Papanikolaou, and P. Dederichs, Phys. 
Rev. Lett. 85, 1088 (2000). 

8 J. Tomfohr and O. Sankey, Phys. Rev. B 65, 245105 (2002). 

9 J. Tomfohr and O. Sankey J. Chem. Phys. 120, 1542 

(2004) . 

10 G. Fagas, A. Kambili, and M. Elstner, Chem. Phys. Lett. 
389, 268 (2004). 

11 E. Prodan and R. Car, Phys. Rev. B 76, 115102 (2007). 

12 E. Prodan and R. Car, Phys. Rev. B in press (2009). 

13 A. Salomon, D. Cahen, S. Lindsay, J. Tomfohr, V. B. En- 
gelkes, and C. D. Frisbie, Adv. Mater. 15, 1881 (2003). 

14 A. Salomon, T. Boecking, C. K. Chan, F. Amy O. Girshe- 
vitz, D. Cahen, and A. Kahn, Phys. Rev. Lett. 95, 266807 

(2005) . 

15 G. Nesher, A. Vilan, H. Cohen, D. Cahen, F. Amy C. K. 
Chan, J. Hwang, and A. Kahn, J. Phys. Chem. B 110, 
14363 (2006). 

16 O. Seitz, T. Bocking, A. Salomon, J. J. Gooding, and 
D. Cahen, Langmuir 22, 6915 (2006). 

17 O. Seitz, A. Vilan, H. Cohen, C. Chan, J. Hwang, A. Kahn, 
and D. Cahen, J. Am. Chem. Soc. 129, 7494 (2007). 

18 A. Salomon, T. Boecking, O. Seitz, T. Markus, F. Amy 

C. Chan, W. Zhao, D. Cahen, and A. Kahn, Adv. Mater. 
19, 445 (2007). 

19 F. Thieblemont, O. Seitz, A. Vilan, H. Cohen, E. Salomon, 
A. Kahn, and D. Cahen, Adv. Mater. 20, 3931 (2008). 

20 D. K. Aswal, S. Lenfant, D. Guerin, J. V. Yakhmi, and 

D. Vuillaume, Analytica Chimica Acta 568, 84 (2006). 

21 C. P. Wade and C. E. D. Chidsey, Appl. Phys. Lett. 71, 
1679 (1997). 

22 P. Allongue, CH.de Villeneuve, S. Morin, R. Boukher- 
roub, and D. D. M. Wayner, Electrochim. Acta 45, 4591 
(2000). 

23 A. B. Sieval, R. Linke, H. Zuilhof, and E. J. R. Sudholter, 
Adv. Mater. 12, 1457 (2000). 

24 A. B. Sieval, V. Vleeming, H. Zuilhof, and E. J. R. Sudho, 



Langmuir 15, 8288 (1999). 

25 A. B. Sieval, A. L. Demirel, J. W. M. Nissink, M. R. Lin- 
ford, J. H. van der Maas, W. H. de Jeu, H. Zuilhof, and 

E. J. R. Sudho, Langmuir 14, 1759 (1998). 

26 L. Segev, A. Salomon, A. Natan, D. Cahen, L. Kronik, 

F. Amy C. K. Chan, and A. Kahn, Phys. Rev. B 74, 
165323 (2006). 

27 A. B. Sieval, B. van den Hout, H. Zuilhof, and E. H. R. 
Sudholter, Langmuir 16, 2987 (2000). 

28 A. B. Sieval, B. van den Hout, H. Zuilhof, and E. J. R. 
Sudho, Langmuir 17, 2172 (2001). 

29 G. Vignale and W. Kohn, Phys. Rev. Lett. 77, 2037 (1996). 

30 G. Vignale, C. Ullrich, and S. Conti, Phys. Rev. Lett. 79, 
4878 (1997). 

31 D. Fisher and P. Lee, Phys. Rev. B 23, 6851 (1981). 

32 H. Baranger and A. Stone, Phys. Rev. B 40, 8169 (1989). 

33 E. Prodan, Phys. Rev. B 73, 035128 (2006). 

34 E. Prodan and R. Car, Nano Lett 8, 1771 (2008). 

35 J. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

36 D. Ceperley and B. Alder, Phys. Rev. Lett. 45, 566 (1980). 

37 N. Troullier and J. Martins, Phys. Rev. B 43, 1993 (1991). 

38 D. D. D. Ma, C. S. Lee, F. C. K. Au, S. Y. Tong, and S. T. 
Lee, Science 299, 1874 (2003). 

39 J. D. Holmes, K. P. Johnston, R. C. Doty and B. A. Ko- 
rgel, Science 287, 1471 (2000). 

40 Y. Zhao and B. I. Yakobson, Phys. Rev. Lett. 91, 035501 
(2003). 

41 X. Zhao, C. M. Wei, L. Yang, and M. Y. Chou, Phys. Rev. 
Lett. 92, 236805 (2004). 

42 H. Scheel, S. Reich, and C. Thomsen, Phys. Stat. Sol. 242, 
2474 (2005). 

43 T. Markussen, R. Rurali, M. Brandbyge, and A. Jauho, 
Phys. Rev. B 74, 245313 (2006). 

44 R. Rurali, B. Aradi, T. Frauenheim, and A. Gali, Phys. 
Rev. B 76, 113303 (2007). 

45 F. Picaud, A. Smogunov, A. D. Corso, and E. Tosatti, J. 
Phys.: Condens. Matter 15, 3731 (2003). 

46 M. Miao, P. van Camp, and V. van Doren, Phys. Rev. B 
54, 10430 (1996). 

47 M. J. Verstraete, P. Bokes, and R. W. Godby, J. Chem. 
Phys. 130, 124715 (2009). 

48 L. Venkataraman, J. Klare, I. Tarn, C. Nuckolls, M. Hy- 
bertsen, and M. Steigerwald, Nano Lett. 6, 458 (2006). 

49 A. Nitzan and M. A. Ratner, Science 300, 1384 (2003). 

50 M. Galperin, M. A. Ratner, and A. Nitzan, J. Phys.: Con- 
dens. Matter 19, 103201 (2007). 



